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Abstract -We study the thermal fluctuation induced interactions between two surfaces containing 
Brownian charges which are held at different temperatures. Using a dynamical form of Debye- 
Hiickel theory implemented within the stochastic equation for the density of mobile Brownian 
charges, we derive expressions for the average force between the two surfaces as well as its variance. 
The latter is found both for the normal, of finite mean, as well as the lateral force, of zero mean, 
between the surfaces. 


Introduction. — Recently there has been intense in¬ 
terest in the out of equilibrium Casimir and/or van der 
Waals interactions between objects held at different tem¬ 
peratures [l]-[^. The analysis of such systems, taking into 
account both quantum and finite temperature effects, is 
carried out using the stoehastie eleetrodynamies formal¬ 
ism of Rytov . This is a very general formalism based 
on Maxwell’s equations in the presence of random cur¬ 
rents which satisfy the quantum fluctuation dissipation 
theorem. The correlation function of currents in the same 
body, assumed to be held at fixed temperature, obeys a 
quantum fluctuation dissipation theorem with the temper¬ 
ature of the body in question. The electrical properties of 
the media in the Rytov formalism are encoded in their 
dielectric functions. 

Here we will study a classical problem from a more mi¬ 
croscopic, and statistical mechanics, point of view. The 
basic model consists of two plates in which a single species 
of charge carrier, with a uniform electroneutralizing back¬ 
ground charge (a so called jellium model), is described 
by an over-damped Langevin equation with an applied 
force generated by a direct Coulomb interaction with the 
charges in the same and opposing plates. In this formal¬ 
ism the dynamical temperature of each plate is well de¬ 
fined and appears as the ratio of the Langevin noise (dif¬ 
fusion constant) to the mobility relating the drift to the 


applied force (the local Einstein relation). Within this 
approximation we compute the average thermal Casimir 
force between the two plates and are able to give explicit 
formulas for the near and far field forces in terms of the 
microscopic parameters of the problem. The far field re¬ 
sults are found to have the same dependence on the plate 
separation L as results derived within the Rytov formal¬ 
ism [T||^ . In addition the model here allows us to compute 
the variance of both the normal and lateral forces in this 
out of equilibrium context and to our knowledge these are 
the first results on force fluctuation statistics for the non¬ 
equilibrium thermal Casimir effect. 

The model. — We will consider a system of two par¬ 
allel planes of area S separated by a distance L as shown 
in fig. §■ Each plane a contains a uniform background 
charge density (per unit area) along with a species 
of mobile ion a, of charge Qa and diffusion constant Dq,, 
obeying the Langevin equation 

= PaDaQa^ + ( 1 ) 

Above is zero mean Gaussian white noise with corre¬ 
lation function {r]ai{t)Vaj{t')) = — t') and E is the 

local electric field in the plane which is generated by the 
electric charge distributions in both planes. The temper¬ 
ature of the plate a is Tc^ = 1//3 q, and thus the electrical 
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plate 2 
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X 

Fig. 1: Schematic of the model, two plates 1 and 2 containing 
Brownian charges of charge qi and q 2 (not shown) respectively. 
The separation between the plates is L and the in plane {x^ y) 
coordinate axes is shown for plate 1. All charges interact via 
the direct non-retarded three dimensional Coulomb interaction. 



mobility of the species a is given by = PaDa from the 
Einstein relation. The statics of the above model within 
the Debye-Hiickel approximation has been studied by a 
number of authors to understand how charge fluctuations 
in neutral colloidal systems give rise to attractive interac¬ 
tions 10 as well as to better understand the high tem¬ 


perature limit of the Casimir force between conductors 
|lll|12]. More recently 13 a stochastic density functional 


approach was applied to the same system in order to see 
how the thermal Casimir force evolves in time between 
two plates having initially no charge fluctuations. Inter¬ 
estingly the relaxation of the thermal Casimir force to its 
equilibrium value turns out to be very different from that 
for a model dielectric based on polarisable dipole dynam¬ 
ics 14 , even if both models can have identical equilibrium 


Casimir forces in the appropriate limits. 

Within each plate a the in plane component of the elec¬ 
tric field can be expressed in terms of the density operator 
of each charged species defined by 


N 


Pa(x) = ^(5(x- Xaj) 


( 2 ) 


i=l 


and is given in the plane a/P hy 


E „/0 = -- [ dx! [V||G(x-x',0)g„//?Pa//3(x') 

^ Js 

+V||G(x - x', L)g^/„p^/„(x')] , (3) 

where V|| denotes the gradient operator in the {x^y) 
planes containing the charges. Here e is the dielectric con¬ 
stant of the system which is assumed to be constant and 
G is the Green’s function obeying 

[Vj^l + ^]G(x -x!,z- z') = -5{x - x!)5{z - z'). (4) 

We also note that in Eq. (§ the uniform background 
charge in both planes do not contribute to the in plane 
electric field. It is useful to define the fluctuations of the 


overall density (per unit area) of species a with respect to 
the mean value via 


Pa(x) = /9„ + n„(x). 


(5) 


In terms of the field nQ,(x) one can write the force acting 
perpendicular to the plates as 

f_^ = JJ dxdx'ni{x)n 2 {x')-^G{x-x',L). (6) 

The i — th component force parallel to the plates, which 
will be on average zero, is given by 



II 


dxdx' ni(x)n2(x')ViG(x 


-x',L). (7) 


The density fields pa in each plate can be shown to evolve 
according to the stochastic partial differential equation 



■ [V||/0a(x,t) -/3a9aEc«/9c«(x,t)] 

+ V| I • [y2Dapa{^,t)Va (x, t)]. (8) 


In the above, the noise t) is a spatio-temporal white 

noise vector field of mean zero and with correlation func¬ 
tion 


(7?„i(x, t)r]pj{x', t)) = dapSijS{t - t')S{x - x'). (9) 

Dynamical Debye-Hiickel Approximation. — As 

it stands eq. ([^ is a non-linear (as the electric field de¬ 
pends on the density) partial differential equation with a 
noise term that depends on the local density field and is 
not amenable to an exact mathematical analysis. How¬ 
ever, one can linearise the equation about the mean value 
of the density to find a dynamical theory which yields the 
Debye-Hiickel high temperature/weak coupling theory as 
its static limit [^. The equation is expanded to first or¬ 
der in Ua in the deterministic part and to zeroth order 
in the noise part (the first order correction to the noise 
is on average zero so this is compatible with expanding 
the deterministic part to first order). The resulting linear 
equation in Eourier (Q) space reads, using the Einstein 
summation convention on repeated indices, 

dtha{Q,t) = -M(Q)a/3%(Q,f) +^a(Q,f), (10) 

where the noise term Ca{Q,t) has correlation function 

= 2(27r)2(5(P + Q)S{t - t')Ra0{Q), 

( 11 ) 

where 

Radio,) = SadDaPaQ^ (12) 

and 

M(Q) = 

q 2 fDi + qiDiPipiGiO-0) 9i92-Di^iPiG(Q, I/) \ 

^ V 9i92il2^2P2G(Q, L) D2 + q^D2p2P2GiO, 0)J ’ 

(13) 
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where G(Q, z) = I2Q is the Fourier transform of G 
taken in the (x, y) plane. From this one can show 17 that 
the steady state correlation function defined by 


{na{Q)M'P)) = (27r)"<5(Q + P)C„/3(Q) (14) 


obeys 

Ma^{Q)C^0{Cl) + Ca^{Q)M^p{Q) = 2i?„^(Q). (15) 


One now has to solve the above linear equation with 3 vari¬ 
ables to obtain the corresponding correlation functions. 
Note that the same formalism could be applied to sys¬ 
tems having more than one charge type; the correspond¬ 
ing equations, while remaining linear, will increase sub¬ 
stantially in complexity. The solution for the steady state 
correlation functions is 


see 


eqs. and @ 


where rUa = PaPa^a/^ inverse Debye lengths (for each 
system confined to a plane). 

The average value of the perpendicular force. — 

From the steady state correlation functions it is now easy 
to compute the average thermal Casimir force between 
the plates in the non-equilibrium state considered here, 
we find that 


{ML)) = 

_ f dQ _ QTiTzmimae-^Q^ 

J (27r)2 (nil + 2Q)(to2 + 2Q) - TOiTO2e“2Qi 

/SiDiTfli + /?2-D2?Tl2 “I" 2(/3iZli + j32D2)Q 

DiTtli + i?2^2 + 2(i?i + D 2 )Q 

( 20 ) 


If one sets Ti = T 2 = T in the above, the second line of 
the above (which is part of the integrand) is equal to P = 
1/T and we immediately recover the equilibrium result 
found in Also if the two systems are at 
different temperatures but are materially identieal in that 
Di = D 2 and mi = m 2 we see that 


(/^(L,Ti,T2))) = 




( 21 ) 


this is a relation first noticed in for planar geometries 
and subsequently shown to be general within the Rytov 
formulation between identical dielectric bodies held at dif¬ 
ferent temperatures [^. 

Of particular interest in this model is the far field limit, 
defined by L ^ where in the equilibrium case 

one recovers the universal thermal Casimir force between 
ideal conductors 13 . For large L the integral in 


eq. (20) is dominated by small Q and one finds 


(/±W) 


SC^{p) ^ T2Dimi-\-TiD2m2 
SttL^ Dimi + D2m2 


( 22 ) 


where (^{z) denotes the Riemann zeta function. When 
Ti = T 2 = T the above becomes independent of and 
ma and we find a long distance force which is universal and 
independent of the microscopic details of the two planes. 
However once the temperatures are different, even the long 
distance force depends on the variables Va = Dama which 
have the dimensions of velocity. The form of eq. (22) is re¬ 
markably simple, but its physical interpretation is far from 
being obvious. We also notice that if we take vi ^ V 2 
by taking a very large diffusion constant Di or large qi 
or Pi the force becomes that of a system with a single 
temperature T 2 at equilibrium. In fact the plane 1 be¬ 
haves like a perfect conductor (due to its ability to react 
instantaneously to an applied electric field) and the ther¬ 
mal Casimir force is independent of its temperature. The 
general behavior of eq. (22) is of the same form found 
for the thermal component of the Casimir force in [TJj^ 
where a dielectric system was considered. Note that in 
this classical model there are no constant bulk contribu¬ 
tions to the Casimir force. This is because here the only 
electrical fields are generated by the charges composing 
the two planes in the systems whereas in a full quantum 
calculation there are additional contributions to the elec¬ 
tromagnetic field coming from blackbody radiation. 


In the near field limit defined by L ^ mp the 

dominant contribution in eq. (20) comes from modes with 
wave-vectors of amplitude Q ^ L~^ and here we obtain 


(ML)) « -5 


/ T2D1 + T1D2 A 

V D 1 + D 2 ) 


mim 2 

IGttL 


(23) 


In the case where Ti = T 2 = T we recover the equilib¬ 
rium near-held result. However we see here that the way 
in which the overall force is weighted with the two tem¬ 
peratures depends on the diffusion constants in con¬ 
trast to the far held limit where it is determined by the 
velocities The results for the average force in the sys¬ 
tem with two plates at diherent temperature show that 
the force in both the near and far held regimes keep the 
same functional dependence on L and that the extent of 
both regimes is still determined by the relative values of L 
and m~^. The weighting of the contribution of each tem¬ 
perature of the force however depends on the dynamical 
variables and We note that the full temperature 
dependence for this system is fully determined by the de¬ 
pendence of mcx on the temperature (which is known), as 
well as the dependence of on the temperature of the 
plate a. 


Force fluctuations. — As both the perpendicular 
force eq. and lateral force eq. 0 are given in terms 
of the density fluctuations, we can see that fluctuations of 
the densities will induce the fluctuations of the forces. 

Using Wick’s theorem we And that the variance of the 
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Cii(Q) 

Ci2(Q) 

C22(Q) 


2Qpi (32 ~ 1)772177126 + (2Q + 7772)(2(i^i + D2)Q + Dimi + 1 ^ 27772 )) 

(2(L)l + D 2 )Q + Di777i + D 2 m 2 ){{ 2 Q + 777i)(2Q + m 2 ) - 777 i7772e-2Q^) 

^ fo) — ‘^QiQ2PiP 2Q e (^11^1 (2(5 + 777i) +/321^2 (2Q + 7772)) 

6(2(L)i + L)2)Q + 1^17771 + D2m2){{2Q + 777i)(2(5 + m 2 ) - 777177726-2^^) 

2 Qp 2 { J 32( fi 2 l pi — 1)777177726 + {2Q + 777 i)( 2 (Di + D 2 )Q + Dimi + ^^27772)) 

{2{Di + D2 )Q + Di777i + D2m2){{2Q + 777i)(2(5 + m 2 ) - 777177726-2^^) 


(16) 

(17) 

(18) 



Fig. 2: Behavior of the variance of the perpendicular ther¬ 
mal Casimir force in the far field limit, as a function of X = 
(Himi — ZI 27772 )/(Hi7771 +II 27772 ) and Y = (Ti — T 2 )/(Ti +T 2 ). 


perpendicular force is given by 


{ fl)c = S 




/ 


dQ 


Cn(-Q)e22(Q)^^^^^ + 


Ci2(-Q)C2i(Q) 




dL 
\ 2 


dL 


dL 


(24) 


Substituting in the correlation function into this expres¬ 
sion then yields 


UlU=sj- 


T 1 T 2 77717772 6 


dQ 

(27r)2 [{{2Q + mi){2Q + m 2 ) - 777i77726-2<3^)2 
7^2(1 - P2/Pi)mim2e-‘^^^ \ 
{2Q + 777i)Hi -j- {2Q + 7772 ) 7^2 / 
d^i(l - Pi/P2)mim2e~‘^^^ \ 


X 1^2(5 + 7771 - 
X ( 2(5 + 7772 - 


(2(5 + 777i)Hi -f (2(5 + 7772 ) 7^2 , 

777^ 7772 


((2(5 + 777i)(2(5 + 7772) - 777 i77726 -2<3^)2 

^ /T2-D1 (2(5 + 777i) + T1H2 (2(5 + 7772) A 
y Hi(2(5 + 777i) + H2(2(5 + 7772 ) / 


(25) 


The variance of the force fluctuations is thus proportional 
to the surface S of the plates, this is typical thermody- 



Fig. 3: Behavior of the variance of the lateral thermal Casimir 
force in the far field limit, as a function ofX = (Himi — 
H 27772 )/(Hi?77i + H 27772 ) and Y = (Ti — T 2 )/(Ti + T 2 ). 


namic scaling for equilibrium systems and persists in the 
non-equilibrium case when the plate temperatures are un¬ 
equal. The existence of such thermodynamic scaling along 
with the absence of any divergences requiring the introduc¬ 
tion of an ultra-violet cut-off, ensures that the variance 


of the force fluctuations must take the form of eq. (26) 


by dimensional analysis. We note that zero temperature 
studies of the quantum fluctuations of the Casimir force 
between two conductors require the introduction of a cut¬ 
off corresponding to the averaging of the force over a finite 
time 18 - 21 . Studies of the thermal fluctuations of the 


force between dielectrics find the same behavior as found 
here with but with a different prefactor [^. The case 
of massless free fluctuating fields with Dirichlet bound¬ 
ary conditions also has an L (plate separation) dependent 
part of the force fluctuation variance which has the same 
form as eq. ^7\ [^ . In both the quantum case and the 
thermal field case the force on a single (or two infinitely 
separated plates) has a non-zero variance; as mentioned 
above this is due to the fact that in these studies there is 
a contribution to the electric field which is independent of 
the charges in the interacting bodies and is generated by 
quantum 18- 21 or thermal 23 fluctuations of the field 


in the surrounding medium. 


The rather complicated expression eq. (25) can be sim¬ 


plified in the near and far-field limits. In the far field limit 
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it is given by 

502 


(/: 


L/c 


2 t^L^ 


3C(3) 


240 8 




V 2 


60 7 


80 


4 7 


(26) 


The average value of the lateral force is clearly zero, 
whereas its connected fluctuation correlator is given by 

(/ll, J||,.)c = ^ L) 

X [Cn(-Q)C' 22 (Q)G(-Q,L) 

-Ci2 (-Q)C'2 i(Q)G(Q,L)]. (30) 


where 0 = (Ti +T 2 )/ 2 , and the temperature and velocity 
mismatches are defined by T = (Ti — T 2 )/(Ti + T 2 ) and 
X = {vi — V 2 )/{vi + '^ 2 )- We thus see that in the far field 
limit {fl)c/{f±)^ ~ L'^/S. 

The behavior of {f]_)c is plotted in Fig. We see that 
{fj_)c is maximum for X = Y = 0, and decreases to zero 
at X = T = ±1. From the above equation we also see 
that in the equilibrium case Ti = T 2 we have 


The calculation of the fluctuations of the lateral force can 
almost be read of from that for the perpendicular force, as 
9G(Q, L)/dL = —QG{Q, L), and one finds in the far-held 
regime that 


(/||,i/||,j)c 


502 

47rL4 


240 8 


-y2 +- X^Y 


8 240 


U^J. 

(31) 


ifY 


2 


) 


c 



(27) 


The term proportional to XT above is negative, thus 
if X and T have the same sign huctuations are sup¬ 
pressed. If aa denotes the ionic radius of the ions in 
plate a and r]'^ the background viscosity, then we see that 
Va = 'p^q^/QeiTiri'^aa by applying the Stokes formula for 
the diffusion constant (assuming hard sphere ions). If the 
viscosity depends only weakly on the temperature then 
so does X. In this regime changing the relative temper¬ 
ature in each plate then either increases or decreases the 
huctuations. 

We have also seen that in the far held limit, when vi = 
V 2 the average equilibrium force is the arithmetic mean of 
two equilibrium systems, one with temperature Ti and the 
other with temperature T 2 . In this case X = 0 and we see 
that 



We thus see that in this case, and in particular the case 
of materially identical plates but with diherent tempera¬ 
tures, the non-equilibrium system has a force which huc- 
tuates less than the corresponding equilibrium one where 
each plate has the temperature 0 = (Ti + T 2 )/ 2 . 

In the near field limit we find that 


{f±iL?)c 


STiT2mim2 

327rL2 


(29) 


Recalling that one must have L m~^ to be in this 
regime and that ttIq, = we see that the fluctua¬ 

tions of the force are independent of the two temperatures 
in this regime. We also see that the dominant fluctuations 
of the near-held force are also independent of the diffusion 
constants of the charges in the two plates. In the near 
held limit we hnd that {f]_)c/{f±)‘^ ^ \ISm\m 2 and so 
the relative huctuations of the force with respect to its 
average value become independent of L. 


In principle there can be correlations between the lateral 
and transverse forces, but in the present case it can be 
verihed there is no correlation. The i = j contribution 
gives the variance of the lateral force in an arbitrary di¬ 
rection. The behavior is plotted in hg. We see that 
(/|^)c is maximum for the plates in thermal equilibrium 
(T = 0), and decreases to zero at |X| = \Y\ = 1. In the 
equilibrium case we hnd from eq. © that 


(e,)2 


(32) 


As is the case for the perpendicular force huctuations, 
when X = 0 we see that the ehect of a temperature diher- 
ence is to decrease the huctuations with respect to that of 
an equilibrium system with the same average temperature. 
An interesting diherence between the lateral force huctu¬ 
ations and the perpendicular force huctuations is that the 
latter depends on the relative signs of X and T through 
its dependence on XT while the former does not. Finally 
the dominant near held contribution to the lateral force 
huctuations is identical to that for the huctuations of the 
perpendicular force given in eq. (29). 

In the case where plate 1 is taken to be a perfect con¬ 
ductor, and so X = 1, we see that in the far held limit 


(/||^)c 


ST 2 /3C(3)T2 
27rL4 V 4 
5TiT2^3 
480L4 


^"(Ti-T2)\ 
240 7 


(33a) 

(33b) 


we see here that while the average perpendicular force in 
this case in independent of the temperature Ti of the per¬ 
fect conductor, its variance and that of the lateral force do 
dependent explicitly on the temperature of the perfectly 
conducting plate. 

Conclusions. — We have studied a simple model 
of Brownian conductor plates giving rise to a thermal 
Casimir interaction between the said plates. The model 
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enables us to study the out of equilibrium behavior of 
thermal Casimir interaction between the plates when both 
plates have different temperatures. We are able to obtain 
explicit expressions for the non-equilibrium attraction be¬ 
tween the plates as a function of the microscopic param¬ 
eters of the model. In particular the far field interaction 
shows how the universal Casimir interaction between con¬ 
ductors is modified by a temperature difference. In con¬ 
trast to the equilibrium case, the non-equilibrium far field 
force depends explicitly on the dynamics of the charges 
through their respective diffusion constants as well as the 
static Debye lengths of both sets of charges. A similar 
dependence is also seen in the near-held force, however in 
both cases the dependence of the force on L is the same as 
that found in equilibrium {l/L^ for the far-held and 1/L 
for the near-held). The non-equilibrium modihcation can 
thus be thought of as giving a diherent non-equilibrium 
Hamaker constant in both regimes. The formalism em¬ 
ployed has also allowed the computation of the huctua- 
tions of both the perpendicular and lateral forces between 
the plates. 

Again the non-equilibrium system has force huctuations 
which have the same dependence on L as the equilib¬ 
rium force, only the temperature dependent prefactors are 
changed. We have also seen that when both systems are 
identical^ in that they have the same Debye lengths and 
diffusion constants (but note that the Debye length de¬ 
pends on the temperature of the plate), the overall force 
is given by the arithmetic average of two equilibrium sys¬ 
tems at both temperatures (as is generally the case for two 
identical bodies in the Rytov formalism !)• However, this 
does not hold for the fluctuations of the force (both per¬ 
pendicular and lateral), which are smaller between identi¬ 
cal plates with different temperatures than the equilibrium 
case where both have the same temperature (Ti + T 2 )/ 2 . 

An obvious extension of this work would be to address 
the effects of inertia on the ion dynamics and also to treat 
the electromagnetic interactions relativistically. Using a 
stochastic density functional theory exists for partially 
damped dynamics, involving both the density fields and 
currents of both species 24 , may render a relativistic cal¬ 


culation feasible. However, the far field results here should 
be relevant up to some length scale, as the near field regime 
is determined by the parameters m~^ which are typically 
nanometric for electrolyte systems. The size of the rela¬ 
tive relativistic correction to the Coulomb potential gen¬ 
erated electric field by a moving charge in one plane in 
the other can be derived from a low velocity expansion of 
the Lienard-Wiechert potential. This correction the has 
order of magnitude R ^ Lvjc? ^ where v is the particle 
velocity and c the speed of light. However for a Brown¬ 
ian particle, both the velocity v and acceleration v have 
correlation functions which diverge at equal times. There 
must however be a microscopic time scale r below over 
which the velocity r is correlated (the ballistic regime). 
Using this, to regularize the calculation of the variance 
of R, gives R ^ LyjDj^j(?. Applying the Stokes Ein¬ 


stein for the ionic diffusion constant D and the estimates 
r = a^pf /t] [^, where a is the ionic radius and p/ is the 
volumic mass density, then suggests that retardation ef¬ 
fects are negligible up to length scales of the order of a 
meter for nano metric ions in water. 
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